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ABSTRACT 


The nature and merits of Tegulometric Frequency Analysis 
(TFA) in determining the presence of weak signals in a noise 
background have been reported in previous works /1-6/7. TFA 
is based on the statistical properties of the signal plus 
noise. Frequency is determined by a zero crossing technique. 
Due to the complex nature of the beating phenomenon in narrow- 
band filtered signal plus noise waveforms, it has not been 
possible to satisfactorily describe TFA in analytical terms. 
This thesis presents the results of an experimental investi- 
gation of the feasibility of a new technique for detection of 
a Sinusoidal signal in noise by defining acceptance parameters 
for amplitude and constancy of frequency. The selection af 
these statistical parameters has been demonstrated to enhance 


the detection of a weak signal in noise and produce an 


effective “processing gain". 
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I. INTRODUCTION 


A. BACKGROUND 

Tegulometric Frequency Analysis (TFA) was developed by 
Dr. George Marmont in 1973 at the Naval Postgraduate School, 
Monterey, California. The analysis determines the presence 
of weak signals in a noise background. The method of analysis 
is based on the statistical properties of the signal plus 
noise, In this method the time between negative to positive 
zero crossings is converted to frequency. A histogram is 
produced which contains the count of the occurrence of each 
frequency as determined by the algorithm used. The theses 
of S. E. Dollar / 1_/ and W. E. Stockslager / 2_/ collectively 
defined the tegule, provided limited mathematical descriptive 
background, and indicated the initial hardware and software 
implementations associated with TFA. A detailed description 
of TFA is contained in two interim reports by Dr. Marmont 
Zz -47, In addition, the processing steps were verified in 
1976 by K. A. Tobin / 5 7. 

The “waxing and warming" nature of a narrow band filtered 
time domain waveform is shown in Figure l. This structure 
was defined as superimposed tegules, The low amplitude por- 
tions of the tegules were observed to randomly contain either 
phase reversals or low amplitude shapes. This fact gave 


support to the heuristic explanation that the observed large 


number of counts associated with weak signals resulted pri- 
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marily from counts occurring during tegule depressions due 


to noise beating against noise. 


B. OBJECTIVE 

The objective of this thesis was to investigate the 
feasibility of comparing processed output data to selective 
amplitude and constancy of frequency acceptance parameters 
to enhance by selection, the detection of a weak signal ina 
noise background. To accomplish this objective the factors 
which influence a relatively large histogram count for a 
signal frequency compared to the histogram counts associated 
with the frequency components of noise had to be identified. 
The major work of this thesis was directed towards identifying 


these factors. 
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II. OVERVIEW OF TEGULOMETRIC FREQUENCY ANALYSIS 


A. HARDWARE 

The particular method used to implement TFA depends upon 
available hardware. 11 data was converted to digital form 
by a twelve bit analog to digital (A/D) converter under the 
control of a minicomputer, The minicomputer is a Digital 
Equipment Corporation (DEC) PDP 11/40. A Time/Data FPE4 fast 
Fourier transform (FFT) microprocessor is attached to the 
PDP 11/40. The TFA algorithm uses the results of the FFT to 
accomplish ideal digital filtering and frequency translation. 
Frequency determination is accomplished in the time domain by 
a zero crossing technique. One of the characteristics of the 
Time/Data PFE4 FET is that it is capable of performing com- 
plex to complex transforms faster than real to complex trans- 
forms. To take advantage of this time saving method one has 
to perform an initial real to complex DFT followed by a 
zeroing of the high frequencies in the upper half of the 
block and a complex to real IFT. Reference 5 contains a more 
detailed explanation of this procedure. Two DEC RKOS5 disk 


drives are used to store the digital data and program routines. 


B. HOW TEGULOMETRIC ANALYSIS WORKS AND WHY 

The basic program to accomplish TFA is labeled HISCAN, 
This program contains subroutines which digitally perform the 
following functions: 


1, narrowband filtering of the time domain waveform. 


at 





*a 


2. detecting positive going zero crossings. 

3. converting the time between positive going zero 
crossings to frequency. 

4. incrementing the count for a detected frequency. 

5. equalizing the counts and thus establishing a histo- 
gram. 

6. limited smoothing of the histogram, if desired. 

Equalization is needed, since for a given time period, 
the counts associated with high frequencies will be larger 
than those of low frequencies. Equalization is accomplished 
by multiplying the frequency dependent word in the histogram 
by a scaling factor whose value decreases with increasing 
frequency. In this way the amplitude of the histogram be- 
comes a measure of the percentage of time that any given 
frequency is observed during the analysis period. 

The histogram is smoothed by a program which adds 
successive and adjacent histogram words to produce a running 
average. The number of frequencies included in this running 
average is entered as an input variable. Figure 2 shows the 
output format of the histogram, The horizontal axis is fre- 
quency and the vertical axis is proportional to the counts 
of the individual frequencies, 

In dealing with digital information, it is convenient to 
establish data blocks in multiples of the number 1024=27°, 
This is referred to as a 1K block. The blocks are labeled BO, 


Bl, B2, etc. Collectively all of the blocks define the memory 
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map for the primary program HISCAN. The data flow in HISCAN 
is shown in Figure 3, which initially appeared in / 5 _/. Here 
a 2K frame of A/D converted data is read from the disk and 
loaded into the 2K block BO. A copy of the data in BO is 
moved to B4. The next sequential 2K frame of data is read 
from the disk and again stored in BO. A copy of the second 
2K frame of data in BO is loaded into the 2K block B9. In 
this manner block BO is used to update the data flow. B4 and 
B9 collectively form the 4K block B8. Processing starts at 
this point. As discussed in the previous section, the first 
step is to get the data into a form where complex to complex 
transforms can be utilized. This is accomplished by per- 
forming a real to complex DFT on block B&. The results are 
stored in block B8 as shown in Figure 3. The upper half of 
block B8, which is block B9, is zeroed and a complex to real 
IFT is performed with the results again being stored in block 
B8. It should be noted at this time that by zeroing block 
B9, the analysis is performed on the lower half of the fre- 
quency components of the data. If an analysis of the upper 
half of the frequency components of the data is desired all 
that is required is a transfer of the contents of block B9 
to block B4 prior to zeroing block B9. If this is done, then 
an analysis of the upper half of the frequency components of 
the data can continue. 

Having performed the above steps, all future transforms 


are processed as complex to complex. This results ina 63.7 
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millisecond savings in time per transform performed. The 
first complex to complex transform is a DFT of block B8 with 
the resultant frequency components again being stored in 
block B8. Repetitive copies of the frequency components of 
the data are moved to block B10 where an ideal digital 
narrowband filter (B2, B15, Bll) is created, This is accom- 
plished by zeroing the frequency components outside of the 
bandwidth of the filter. A complex to complex IFT is then 
performed on block B10. The 4K data block in B10 now re- 
presents a narrowband filtered time signal and, if plotted, 
displays the waxing and wanning nature shown in Figure l. 

- To reduce significantly the end effects associated with 
the discrete inverse fourier transfer, block B5 is defined 

i as the center 2K data points of the 4K block B10. The upper 
1K data points and the lower 1K data points are discarded. 
A rather ingenious method of data flow will recover the in- 
formation within these points. This will be discussed in the 
following paragraph. The processing continues from here by 
applying the zero crossing detection routine, called TEGULO, 


to block B5. TEGULO determines the instantaneous frequency 


between positive going zero crossings. It compares this 

-_ calculated frequency to the known frequencies contained with- 
in the narrowband filter. If the calculated frequency is 

“ within the band, its associated count stored at an address 


in block B7 is incremented, Notice that block B7 is 4K in 


length. 
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To evaluate the frequency components from the minimum to 
the maximum, the narrowband filter is made to scan across the 
band of interest. Accordingly, after TEGULO has been run, 
HISCAN loops back to move a copy of the frequency components 
from block B8 to block B10. The narrowband filter is rede- 
fined so that the lowest frequency element is dropped from 
the filter and the next higher frequency element is added to 
the filter. This process repeats itself until all frequencies 
of interest have been scanned by the narrowband filter. HIS- 
CAN then loops back to load the next 2K frame of sequential 
data from the disk. Recall that TEGULO operates on the 
center 2K data points of the narrowband filtered time signal 
and discards the lower and upper 1K data points. The informa- 
tion contained within these 1K data points is processed as 
follows. Block BO still contains the data loaded into Block 
B9 from the previous loop. By loading a copy of block BO 
into block B4 this information is processed on this pass of 
the loop. Block BO is of course loaded with the next se- 
quential 2K data points from the disk and a copy of this is 


moved to block B9 before processing begins. The sequence is 


shown below. 








Processed Data Blocks 
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* Numbers refer to sequential 1K segments of data. 


Observe that the processed data blocks are in time sequential 





order. This process continues until the specified number of 
2K data frames have been read from the disk. When all the 
data has been read from the disk, a routine called EGUALI is 
applied to the histogram block B7 to equalize the counts for 
high and low frequencies. After EGUALI the histogram is 
smoothed by a routine called RUNAVG. The histogram is now 
ready to be plotted and appears as previously shown in Figure 
oe | 

The above is an outline of the mechanics of TFA. Recent 
studies have identified certain key items which have a marked 
influence on the performance of TFA, These items are discussed 
below and will be presented in more detail in a future report 
to be published by Dr. Marmont. They are outlined here with 


his permission, . 
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To improve the processing capability of TFA a method of 
reducing the noise power competing with the signal power was 
developed. Initial discussion of this effective noise reduc- 
tion method appeared in reference 5 and subsequent improvements 
in this method were reported in reference 6. The method con- 
sists of resampling the data. The data to be analyzed is 
sampled at a rate of 512 Hz. The samples are stored ina 
block 4K in length and therefore represent 8 sec in time. 

This is shown in Figure 4. When the DFT is performed the fre- 
quency resolution is 1/8 Hz per DFT frequency element. The 
DFT frequency element consists of two words, the first word | 
is the real component and the second word is the imaginary 
component of the Fourier series coefficients. The DFT block 
represents frequencies from 0-256 Hz. The data will be re- 
sampled by a software routine called RESAMP. This routine 
digitally resamples the data at a rate of 64 Hz or 1/8th of 
the original 512 Hz sampling rate. Prior to resampling, how- 
ever, the Nyquist criterion must be satisfied. To accomplish 
this, the data must be digitally filtered to a 32 Hz band. 
This 32 Hz band may be located anywhere within the DFT block. 
For example consider the case where the 32 Hz band from 48 Hz 
to 80 Hz is to be analyzed. 48 Hz is represented by the 

384th DFT element and 80 Hz is represented by the 640 DFT 
element. The elements from 384-640 are moved down in the DFT 
block such that element 334 becomes the first element above 
the DC element. In the process of basebanding the frequencies 


of interest, consideration must be given to whether the 
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frequency element is odd. This is discussed in the following 
paragraph. The DC element and the upper 7/8th of the DFT 
block are then zeroed. An IFT of the block is then performed 
resulting in 4K time samples representing 8 seconds in time 
and containing frequencies from 48-80 Hz. This is shown in 
Figure 5. In order for the linear interpolation to be valid, 
there must be at least 10 sample points per cycle (assuming 

a sinusodial waveform). Recall that on resampled data a 32 
Hz band of frequency is basebanded. The 32 Hz band, however, 
satisfies the Nyquist rate giving only 2 sample points per 
cycle at the high end of the band. To obtain 10 sample 
points per cycle the entire rangeof the scan filter must be 
no greater than zero to 6.4 Hz (32/5). As a comparison, the 
bandwidth of the narrowband filter is on the order of 10 fre- 
quency elements. For resampled data this is 10/64 Hz. Ob- 
serve also in Figure 5 that the fractional word count between 
zero crossings is obtained. In order to maintain with inter- 
ger arithmatic the numbers are multiplied by 2" during the 
calculation process. Knowing the time between the words, the 
period and, hence, the frequency can be calculated. The 
frequency is then represented as an address in the histogram 
block. This method provides a 10 fold improvement in the 
frequency resolution. For example, consider the case where 
the data is resampled at 64 Hz. For the 4K data block the 
conventional resolution is 64 DFT elements/Hz, In TFA, how- 


ever, the requirement of 10 sample points per cycle sets the 
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highest frequency at 6.4 Hz. The 6.4 Hz band represented by 


the histogram is stored in a 4K block giving 640 histogram 





points/Hz. 

No discussion of TFA would be complete without giving 
thought to the construction and interaction of the scan filter. 
The filtering is substantially valid if there are at least 
10 consecutive frequency elements in each rectangular window, 
provided the two 1K end blocks of the resultant 4K IFT block 
are discarded. This substantially reduces the end effects 
inherent with the DFT. The effect of the scanning is to 
average the tendency of the addition of several sinusoids to 
"pull" the zero crossing measured frequencies to the center 
frequency. This gives a very linear frequency scale for the 
histogram block. The filter is made to scan by dropping the 
lowest frequency element and adding the next higher frequency 
element. In resampled data this is a aieaes of 1/64 Hz. This 
method gives equal exposure to all the elements. As HISCAN 
reads the next 2K block of data from the disk, the noise will 
be quite different from the previous 2K block of data in 
magnitudes, frequencies, and phases. Each successive sinusoidal 
signal in the next new block, however will be substantially the 
same, except for phase. Hence the successive scans average 
a mixture of noise components and possible signal in various 
combinations of amplitude and phase with respect to frequency. 
Due to the small incremental change in the scan filter (1/64 


Hz for resampled data) the frequencies of the components 
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change gradually. The range of frequencies for a given scan 
number in the case of resampled data is very small (10/64 Hz). 
Although the incremental change and the number of elements 
included in the scan filter are small, the tegular nature of 
the IFT's for successive incremental scan positions is experi- 
mentally determined to differ markedly. This is shown in 
Figure 6 where two IFT's, differing only by one DFT element 
are plotted. Due to the complex nature of the beating phe- 
nomenon and the averaging performed by the scan filter, a 
mathematical analysis of TFA is impractical. An experimental 
approach to determine the important factors is efficient and 


effective. 
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III. REVISIONS TO TFA 


To investigate the nature of the Tegules and the feasibil- 
ity of a selective comparison for processed output data to 
amplitude and constancy of frequency input parameters to en- 
hance a weak signal the subroutine which detects the zero 
crossings (TEGULO) was modified. Recent unpublished studies 
had shown that the input amplitude and constancy of frequency 
could be combined to define a selective criteria. The ampli- 
tude and constancy of frequency of the associated computed 
zero crossings could be compared to the established criteria 
to determine if the histogram count should be incremented. 

To observe the effect of both high and low amplitudes two new 
versions of TEGULO were written. One version (called HISCAN- 
DEL) required the amplitudes to be above a callable percentage 
of the maximum amplitude to be considered for inclusion in the 
histogram and the other version (called HISCAN-XXX) required the 
amplitudes to be below a callable percentage of the maximum 
amplitude to be considered for inclusion in the histogram. 

Both versions established a running average of the change in 
frequency (Af) corresponding to the periods between successive 
nagative to positive zero crossings. The length of the running 
average could also be changed by a callable parameter. The 
instantaneous Af had to be within the selected range of the 
running average to be considered in the histogram. The pro- 


grams were set up to show graphically which zero crossings 
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were included in the histogram. To reduce the end effects asso- 
ciated with the DFT, only the center 2K samples are resampled. 
The routine RESAMP selects every eight sample and stores it in 
a new 2K block. It takes eight passes to the loop to fill the 
new 2K block. When the storage block is filled it represents 
32 seconds in time, and it is transferred to the disk for 
permanent storage. The data flow in RESAMP is the same as that 
in HISCAN to insure that there is no loss of data. Resampling 
thus performs a time compression. The effect of resampling is 
to increase the frequency resolution and to decrease the noise 
per DFT element. This is apparent when thought is given to 
the data flow HISCAN. Recall that HISCAN established the 
4K block B8 (Figure 3). With resampled data, which has a 
sampling rate of 64 Hz, the 4K block B8 represents 64 seconds 
in time. When DFT is performed the frequency resolution is 
1/64 Hz. Therefore the frequency resolution has been in- 
creased by a factor of 8 and the noise per DFT element has 
been reduced by a factor of 8. This increases the effective 
S/N ratio. The numerical value of this increase is developed 
in the section dealing with Test Procedures and Results. 

In the routine RESAMP when the lowest frequency of interest 
to be baseband converted is odd, the time domain representation 
will be 180 degrees out of phase on every other pass through 


the loop. This problem is solved by negating both the real 
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and the imaginary Fourier coefficients for every odd number 
loading of a 2K word data frame, 

The zero crossing routine, TEGULO, performs a linear inter- 
polation between the negative to positive data words at each 


side of the zero crossing. There are four possible outcomes: 


Within Af Range ‘Within amplitude range Graphical Display 
(CRT and/or Plot) 


No* No Af | x 


No Yes Af | 


marker 1 
£ x 


Yes* No A 


Yes Yes = 


| marker 2 


* "DON'T CARE" condition 
marker 1| < marker 2 


Figure 6 shows a typical result. 
With the above modifications to HISCAN, the influence of am- 


plitude and constancy of frequency boundaries on enhancing 
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? the detectabiliby of a weak signal in a noise background 


could be experimentally determined, 





672 


Starting with element 114 


Starting witn element 115 


Figure 6 


IFT of 10 DFT Elements of Resampled Data 
Frequency Content 49.781-49.938 Hz 
S/N = -16 db (1 Hz band) 








IV. THE DATA 


To investigate the nature of the tegules, test signals 
of various signal to noise ratios were produced. The source 
of the noise was a Gaussian noise generator, The RMS voltage 


value of the noise was observed on a true RMS voltmeter. The 





signal was adjusted to give the desired signal to noise ratio. 
The noise and the signal were combined by an analog adder and 
passed through a 4th order Butterworth filter to limit the 
frequency content from 0-200 Hz. The analog signal plus noise 
waveform was digitized by a twelve bit A/D converter. 

The above method was found to be adequate for large S/N 
ratios. When a low S/N ratio was desired the signal could 
not be adequately represented. Due to the dynamic range 
of the 12 bit A/D converter, a signal at S/N = -20db was 
represented by only three bits. This can readily be seen by 
assuming the maximum noise voltage to be roughly 5 times the 
value of the RMS voltage. With an RMS noise voltage value 
of 0.4 volts in a 200 Hz band the A/D converter can have a 
dynamic range of +2 volts without experiencing any clipping. 
The required voltage level of the signal can be caiculated 


ie: from: 


: o,. = 0.4 volts (RMS value of the noise in a 200 Hz2 band) 


o = 0.16 








——- = oo Reduce Bandwidth to 1 Hz 


0.0008 Watts/Hz (assume 1 ohm resistor) 

For a signal to noise ratios of - 20db 

S = 0.000008 watts 
The amplitude for a sinusoidal signal would then be 0.004 
volts. With a 16 bit computer two volts would be represented 
by the number 32767. Therefore: 


0.004 volts = 66 


10 
= 102, (octal) 
. = 1 000 010, (binary) 


With only 12 bits for the A/D converter, the lower 4 bits are 
lost and the signal is represented by only 3 bits. This is- 
not sufficient and a method of utilizing the full 7 bits to 
represent the signal was sought. 

To achieve the maximum precision, the signal and noise 
were separately converted to digital form at a sampling rate 
of 512 Hz and stored on different disks. A routine called 
SIGMA was used to calculate the RMS value of the noise voltage. 
SIGMA transforms the noise to the frequency domain and sums 


the square of the real and imaginary Fourier coefficients to 


. find the variance of the noise. The algorithm utilized is 


1 1 . 
oc = ey 4 5 om (a +b) 
0 
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total number of sample frames 


x 
oy 
@ 
RR 
@ 
s 

Hl 


N = total number of frequency elements per frame 
a_, b. = Fourier series coefficients 


variance of the noise (zero mean) 


Q 
It 


This software calculation compared favorably with the observed 
meter reading for the RMS value of the noise. 

With the RMS value for the noise known a program called 
TONOTS was used to scale the signal voltage and sum the noise 
and the scaled signal. In this way seven bits of significance 
were realized. This data was then resampled by software 
utilizing the program KESAMP. As daiscussed in section II this 
resulted in an effective sampling rate of 64 Hz. The frequency 
resolution was increased by a factor of 8 to 1 uz and the 


64 
noise per DFT element was reduced by a factor of 8. 
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V. TEST PROCEDURES AND RESULTS 


The test procedures consisted of varying the input para- 
meters to HISCAN for the same data and observing the outputs. 
The modified versions of HISCAN had eight input paramters as 
follows: 

PO = starting DFT element 


Pl = the number of DFT elements in the narrowband 
scan filter 


P2 = the starting disk address 

P3 = the number of 2K reads from the disk 

P4 = the number of points to be used in RUNAVG 
P5 = the amplitude factor 


-* P6 = the number of words to be included in tne 
running average of Af 


P7 = the allowable deviation of the instantaneous 
frequency trom the running average of Af. 


PO through P4 remained constant while P5 through P7 were 
varied. P5 was expressed as a fractional part (1/16th) 

of the maximum amplitude. Depending on which version of 
HISCAN was used, P5 established a level above or below which 
the amplitude associated with a zero crossing had to be in 
order for its frequency to be counted in the histogram. P6 
was expressed as an address offset over which the running 
average Af was estaplished, Recail that block B5 in Figure 3 
contains the narrowband filter time waveform, TEGULO is 
applied to block B5. The frequency associated with two con- 


secutive negative to positive zero crossings is computed. 


33 





The change in frequency Af is computed by subtracting the 


frequency associated with the previous negative to positive 





zero crossing pair from the frequency assoicated with the 
current negative to positive zero crossing pair, The value 
of Af is stored in a new block (B23) in a location relative 
to the location of its associated zero crossing vair in block 
B5. B23 is initially zeroed. In this way, the running 
average can be computed by summing the contents of B23 within 
the range of addresses specified by P6 and dividing by the 
number of zero crossing contained within this range. This is 
done for each negative to positive zero crossing pair. The best 
maximum look back established was equal to 1/4 the block 
length of B5 which was 2K. Expressed as an octal address, 


P6 could vary between 0-4000 A vaiue of 1000, was experi-~ 


3° 
mentaliy determined to give the pest results. ‘This corres- 
ponds to using 1/4 of the block to establish the running 
average of Af. 

P7 defined the allowable instantaneous frequency divia- 
tion (Af) from the running average of Af. P7 was expressed 
as an address offset also, In this case nowever, P7 was 
entered as the allowable offset number of bytes. Here a byte 
is equal to 8 bits of the 16 bit computer word. For example, 
if a value of 10, for P7 was used this reduces to 5 computer 


words, For resampled data having 640 computer words per Hz, 


Af = 1/128 Hz, 
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A starting point was established based on past observa- 
tions. These observations indicated that the version of HIS- 
CAN, called HISCAN:DEL, that retained zero crossings associa- 
ted with amplitudes that were above a callable amplitude level 
produced better results than HISCAN*XXX which retained zero 
crossings associated with amplitudes below a callable ampli- 
tude level, With this starting point, the parameters P5, P6 
and P7 were then varied. Outputs for which there were marked 
differences were examined on a CRT to correlate tne zero 
crossings retained with the tegular nature of the waveform. 
The format for the CRT diaplay is shown in Figure 7. Other 
outputs which appeared interesting for one reason or another, 
i.e. no change, were also studied. 

The results of the tests confirmed the preliminary obser- 
vation that the version of HISCAN (HISCAN*DEL) which include 
only data apove an amplitude ievel in the histogram produced 
better results than tne version of HISCAN (HISCAN:XXX) which 
included only data below an amplitude level in tne histogram. 
HISCAN-DEL also produced better results tnan the original 
version of HISCAN which did not contain an amplitude or con- 
stancy of frequency criteria. A typical, result is shown in 
Figure 8. This result appeared contrary to the past heuristic 
explanation that tne majority of counts for a weak signal in 
a noise back ground occurred during teguie depressions due 
to noise beating against noise. Figure 9 has been choosen 


to demonstrate the relative order of magnitude of the zero 
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Counts not accepted . 
for nistogram due 

to frequencies being 
outside of scan filter 
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counts not accepted count not accepted 
for histogram due for histogram due 
to Af to amplitude 
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counts accepted for 
histogram 


Figure 7 


Observation of Tegules with Selective Amplitude 
and Constancy of Frequency Acceptance Parameters 








LE 


49.063 Hz 


50.000 Hz 


Histogram trom original version of HISCAN 


S 


S/N = -16 db (1 Hz band) Resampled data 
16 2K Reads from disk (Acquisition time: 8.533 min) 
= 50 Hz 
= 


Figure 8 (a) 


Comparison of Histograms For Various Amplitude 
And Constancy of Frequency Input Parameters 





51.000 Hz 








Note: Larger relating number of 
counts associatea with tne 
a Signal (50 Hz) 


8E 


49.063 Hz 50.000 Hz 51.000 Hz 


Histogram from HISCAN-DEL. Frequency counted for 


Zero crossings whose assoicated amplitudes were > 1/2 maximum 
amplitude : 


Data, including acquisition time, is tne same as for 
Figure 8(a). 


Figure 8(b) 
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a 
ue Note: Relative large numper of zero crossing counts to be 
included in the nistogram. 


Figure 9 


Histogram Counts Witnout Selective Comparison 
Of Processed Output Data 


crossing counts included in the histogram when the amplitude 
factor has been set to 0 in HISCAN:‘DEL (all amplitudes are 
considered for inclusion in the histogram). This is analogous 
to earlier versions of HISCAN. By comparison Figure 10 shows 
the relative reduction in the order of magnitude of the zero 
crossing counts included in the histogram when the data is 
selectively compared to the input parameters. If the reduc- 
tion in counts in the nistogram is due to "garbage" being 
discarded, tne signal is ennanced. These are typical results 
and demonstrate the potential for a great processing gain, 
i.e., enhancement of the signal counts, provided tne zero 
crossings not considered for count nave a low probability of 
containing the signal. Tnis propability cannot be determined 
analytically. The computer must be used to obtain the num- 
erical results. Time was not available to compute this pro- 
bability and to include it in this thesis. A proposed method 
for determining this probability 1s includea in the following 
section. 

Evaluation of tne test results revealed that it was possible 
to decrease the overall number of counts in the histogram by 
a comparison of the processed output data to an input ampli- 
tude and constancy of frequency criteria and at the same time 
to enhance the signal, This is shown in Figures 11 and 12. 
In Figure 11 the input parameters were 


P5 = 12 (5/ 8th) P6 = 1000 Pi = 10, 


8 8 7 
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Frequencies counted for zero crossings whose 
associated amplitudes > 5/8th of maximum amplitude 
and whose Af < 1/128 Hz 


TP 


NOTE: Relative reduction in zero crossing counts 
to pe included in the histogram 


Figure 10 


Histogram Counts With Selective 
Comparison Of Processed Output Data 








49.063 Hz 


50.000 Hz 


Histogram from H1SCAN:DEL. 
S/N = -16 db (1 Hz band) Resampied data 


16 2K Reads from tne disk (Acquisition time 8.533 min) 
Signal = 50 Hz 


Figure 11 
Histogram For Amplitudes > 5/8th Maximum 
Af < 1/128 Hz 


51.000 Hz 
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4——_ NOTE: Enhancement In Signal 
, Counts 


49.063 Hz : 50.009 Hz - 51.000 Hz 
1. 





Histogram; from H1iSCAN+ DEL : eB 
‘S/N = -16'db (1 Hz band) Resainpled data 


16 2K Reads from thé disk (Acquisition Ped 8.533 min) 
S = 50 Hz 


Figure 12 


Histogram For Amplitudes > 6/8th Maximum 
Af < 1/128 fiz 








While in Figure 12 the input parameters were 


P5 = 14 (6/8th) P6 = 1000, 


For this particular case an increase in the amplitude Level 
above which the counts were included in tne histogram, im- 
proved tne histogram. The selection process therefore re- 
sulted in a "processing gain". In cases such as this, where 
a high count for a valid signal occurred, examination of the 
zero crossing plots (Figure 7) indicated the following: 

1. Af was small over an appreciable interval. 

2. The amplitude was relatively large. 
The fact that the large amplitudes enhanced the signal was 
confusing at first. The heuristic explanation of how TFA 
works has always included the concept that the signal was more 
likely to be observed during the tegule depressions. In the 
case of resampled data, the effective signal to noise ratio 
per 10 DFT spectral elements used in each scan is increased. 
In the previous section it was shown that for a S/N = -20 db 


(in a 1 Hz band) 


Vota = 0,002828 volts RMS 
S 2 

N _ WATTS 
as = 0.0008 as 


As a result of resampling the frequency resolution is in- 
creased to 1/64 Hz. Since a scan width of 10 DFT spectral 
elements was used, the scan bandwidth is reduced to 10/64 Hz 


from 1 Hz. 
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Therefore (for S/N=-20db before resampling), after resampling: 





4 2 WATTS 10 
on = 0.0008 az x 64 Hz 
2 = 0.000125 
on . 
On = 0.01118 
Verte = 0.002828 RMS 
V 
SIG = 0.2529 
On | 
Ss = -_ 
5 = - 11.94 db 
as 


This is an 8.06 db improvement in the effective S/N ratio. 

te For data originally at S/N of -16db(in a 1 Hz band) which was 
the case for the majority of the test runs, the effective 
signal to noise ratio was increased to - 7.94 db. 


Therefore 
= 0.401 (effective) 


rather than the voltage ratio of 0.158 associated with ~-16 db. 
Reference 6 reported on the extremely short acquisition 
- time required for TFA as compared to systems presently used 
by the Navy. This short acquisition time is one of the sa- 
lient features of TFA, Its importance cannot be overemphasized. 


To demonstrate this capability, the acquisition time for 


Figures 8, 11 and 12 (8.533 minutes for -16 db S/N ina l Hz 





band) was reduced by a factor of 1/2 to 4.267 minutes and the 
analysis repeated. The results are shown in Figure 13. The 

50 Hz signal still received the largest number of counts in 

the histograms of the HISCAN’DEL version. Future work needs to 
be directed towards establishing a minimum acquisition time 

for a given signal to noise ratio. 

In summary the results demonstrated that it is feasible 
to compare the processed output data to amplitude and con- 
stancy of frequency input parameters to enhance the signal. 
For test data at a S/N = ~ 16 db (1 Hz band), the histogram 
counts assoicated with the signal was enhanced when the am- 
plitude was 5/8ths or 6/8ths of the maximum amplitude and 


the change in frequency was small over an appreciable inter- 


val. 








49.063 Hz : ’ 50.000 Hz 


Histogram from original version of HISCAN > 


S/N = -16 db (1 Hz band). “Resampled data 
8 2K reads from the disk (Acquisition time 4.267 min) 
Signal = 50 Hz . 


Figure 13 (a) 


Reduced Acquisition Time Processing 


51.000 Hz 





49.063 Hz 50,000 Hz 51.000 Hz 


Histogram from HISCAN:DEL. Frequencies counted 
for zero crossings whose amplitudes were > 1/2 maximum 
Data is the same as for Figure 13(a) 


Figure 13(b) 








Histogram from HISCAN. DEL. pracubaeies counted 
for zero-crossings whose amplitudes were-->-' 5/8 maximum 
S/N = -16 db (1 Hz band) Resamplea data . 

8 2K reads from tne disk (acquisition time 4.267 min) 
Signal = 50 Hz 


Figure 13 (c) 
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49.063 Hz 50.000 Hz 51.000 Hz 


OS 


Histogram from HISCAN’DEL., Frequencies counted 
for zero crossings whose amplitudes were > 6/8 maximum 
Same data as Figure 13(c) 
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Figure 13(d) 





VI. CONCLUSIONS AND RECOMMENDATIONS 


The test results demonstrated that it is feasible in TFA 
2 to enhance a weak signal in a noise background by selecting 
input parameters against which the results of processing are 
to be compared for acceptability. In this context then, selec- 
tive comparison of processed output data results in a "pro- 
cessing gain". The fact that better results are obtained 
when amplitudes of 5/8ths or 6/8ths of the maximum are sel- 
ected for the histogram counts can be explained in the fol- 
lowing manner. By resampling the data, it was shown that the 
Signal to noise ratio in the scan band is increased by a 
factor of 8.06 db. Therefore the RMS signal to noise voltage 
_ ratio in the scan band for an original signal to noise ratio 
of -~16 db in a1 Hz band is increased from 0.158 to 0.401. 
The scan band consisted of 10 complex spectral elements. 
Assuming for simplicity a signal with no DFT leakage, each of 
the 10 complex spectral elements represents noise with the 
exception of one element which represents a combination of 


Signal and noise. 





As the narrowband filter scans across the frequencies of 
interest, the noise changes each time the lowest DFT element 
is dropped and the next higher DFT spectral element is added. 

. The signal, however, does not change, The frequency, ampli- 
tude, and phase of the signal are constant for the particular 


read of data from the disk, When a new block of data is read 
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from the disk, the noise is of course different, but the 
signal differs only in phase (assuming a constant source and 
neglecting effects of the transmission medium). Accordingly, 
when the narrowband filter of 10 complex spectral elements 
scans across the signal there are ten scans during which the 
noise is changing but the signal is not changing. Since the 
effective RMS signal to noise voltage ratio in the scan band 
is 0.401 (for -16 db S/N) and since the amplitude of the 
noise is changing, it is highly probable that the signal am- 
plitude will be greater than the noise amplitude for a num- 
ber of scans. Under these conditions, when the signal ampli- 
‘tude is higher than the noise amplitude, the signal will 
"capture" the noise. This is the same capture effect asso- 
ciated with FM transmissions. The high amplitude frequency 
component will dominate. This will result in several se- 
quential negative to positive zero crossings being at or near 
the signal frequency. With this in mind, the desirability of 
acceptable processed data also based on a relatively slow rate 
of change of frequency becomes apparent. The signal frequency 
as measured by zero crossings will receive a relatively larger 
number of counts than the noise components. The concept of 
"Capture" appears to be a reasonable explanation to the 
correlation of certain input parameters with the histogram 
output. 

It is recommended that since Reference 6 has already de- 


monstrated the capability of TFA to detect weak signals in 
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the ocean environment with short acquisition times, that a 
prototype tegulometric frequency analyzer be developed and 
tested. As a follow-on to this thesis, it is recommended 

that the zero crossing routine TEGULO be modified to compute 
an “enhancement index" of a signal as a function of amplitude 
and constancy of frequency. One possible way to compute this 
enhancement index is to record the noise and signal on separate 
disks, as previously discussed, and to resample the two disks 
separately. Before scaling and combining the signal with the 
noise, HISCAN, with selected input parameters, can be run on 
the noise. The histogram would indicate how many counts at 
the signal frequency were due to noise alone. After combining 
the signal and the noise, the histogram with the same input 
parameters would indicate the number of counts at the signal 
frequency due to signal plus noise. The difference between 
these two counts would be a measure of the enhancement index 
of the signal. The probability of a signal being enhanced for 
a particular amplitude criteria and a given Af would require 
computing the enhancement index for the case where all am- 
plitudes are considered for analysis. This would be followed 
by computing the enhancement index for amplitudes above and 
below the amplitude level criterion. The ratio on the en- 
hancement index above the amplitude level to the enhancement 
index for all amplitudes would be an indication of the pro- 
bability of the signal counts being contained in the ampli- 


tudes above the level. In a similar way the probability of 
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the signal counts being contained in the amplitudes below the 
level can be computed. As a graphical aid, the plotting routine 
can be modified to indicate which counts for the histogram are 
within a specified frequency range of the test signal. The 
above modifications would supply numerical data from which 
the optimum selection of the input parameters and a valid ex- 
planation of the workings of TFA could be obtained. 

Work performed on test data with a S/N ratio of -20 db 
(1 Hz band) revealed that the counts associated with the sig- 
nal were not always the maximum count on the histogram. 
Evaluation of sequential processing results highlighted the 
random nature of the noise and the constancy of the signal 
frequency. It is recommended that for extremely low S/N 


ratios a correlation technique be developed for TFA. It is 





felt that a correlation capability would greatly improve the 
sensitivity of TFA to detect a signal at weaker than -20 db 
S/N ratio (1 Hz band). The implementation of this correlation 
feature would not be difficult. One possible way of achieving 
a type of correlation would be to define two blocks (double 
precision) each parallel to the histogram block. One block 
would contain a product and the other block the sum of the 
products, The first histogram block would be moved to the 
product block. The product block would then be multiplied 

by the histogram block for the next section of data with the 
double precision result being stored in the product block. 


The product block would then be normalized and shifted down 
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to a single precision word. The product block would be added 


to the sum 
have to be 
after each 
tions made 


promise in 


block (initially all zeros), Again the sum would 
normalized and shifted down to single precision 
addition to prevent overflow. With the observa- 
to date this correlation like function holds great 


increasing the sensitivity of TFA. 
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